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P ■ ABSTRACT 

(N 

^ I We report on submillimeter and infrared observations of 28 radio-selected galaxies 

in the Spitzer First Look Survey Verification field (FLSV). All of the radio-selected 
galaxies that show evidence for emission at 850/im with SCUBA have Spitzer counter- 
parts at 24^m, while only half of the radio-selected galaxies without 850//m emission 
have detectable counterparts at 24//m. The data show a wide range of infrared colors 
(S70/S24 < 5-30, S8/S3.6 < 0.3-4), indicative of a mixture of infrared- warm AGN and 
^f^ ', cooler starburst dominated sources. The galaxies showing 850//m emission have Spitzer 

^^O I flux densities and flux density ratios consistent with the range of values expected for 

,—1- ' high-redshift (z = 1-4) ultraluminous infrared galaxies. 

o' 

^S^ • Subject headings: galaxies: active — galaxies: evolution — galaxies: formation 

galaxies: starburst 
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Introduction 



^ . The Spitzer Space Telescope provides us with the exciting opportunity to study the high- 

H ! redshift universe at mid&far- infrared wavelengths. The IRAS mission first uncovered the presence 

of infrared luminous galaxies in the local universe (Neugebauer et al. 1984), and the sub-mm/mm 
surveys with SCUBA and MAMBO have highlighted the importance of ultraluminous infrared 
galaxies (ULIRGs, > 10^^ Lq) at high redshift (e.g.. Small, Ivison, & Blain 1997; Hughes et al. 1998; 
Bertoldi et al. 2000). Studies of high-redshift ULIRGs are important for our general understanding 
of galaxy evolution since they are responsible for a significant fraction of the total energy generated 
by all galaxies over the history of the universe (e.g., Blain et al. 2002). 
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The recent spectroscopic studies of the sub-mm galaxy (SMG) population show that the red- 
shift distribution peaks at z ~ 2-3 (Chapman et al. 2003b, 2004), and that the population is 
comprised of starbursts, AGN, and composite AGN+starburst systems (Ivison et al. 1998, 2000; 
Frayer et al. 2003; Knudsen, van der Werf, & Jaffe 2003). Even though many SMGs show the 
presence of AGN, the molecular CO line emission (Frayer et al. 1998, 1999; Neri et al. 2003) and 
X-ray data (Alexander et al. 2003) are consistent with the majority of the infrared emission for the 
population arising from star formation. The Multiband Imaging Photometer for Spitzer (MIPS, 
Rieke et al. 2004) allows us to directly measure the infrared colors and constrain the fraction of 
infrared-warm AGN dominated versus infrared-cool starburst dominated SMGs. 



2. Observations 

Before the launch of Spitzer, we identified potential SMGs in the First Look Survey Verification 
field (FLSV) by selecting radio sources with faint optical counterparts, following previous successful 
selection techniques (e.g., Cowie, Barger, &: Kneib 2002; Chapman et al. 2003a). We used deep 
Westerbork 1.4 GHz radio data (rms= 9;uJy, Morganti et al. 2004) and deep optical NOAO R- 
band data (3a = 26.4 mag, Fadda et al. 2004a) to derive a list of candidate sources for follow-up 
observations with SCUBA. In the Spring of 2003, we observed 28 galaxies at the James Clerk 
Maxwell Telescope (JCMT) using the SCUBA two-bolometer photometry mode, achieving rms 
levels of 2~3 mJy at 850;um. 

The Spitzer observations were taken as part of the Extragalactic First Look Survey (FLS)^. 
The 28 galaxies in the sample are located within the 0.25 degree^ of the FLS verification field 
(FLSV) and were observed with both the InfraRed Array Camera (IRAC, Fazio et al. 2004) and 
MIPS. The data presented here have effective integration times of 480 s, 336 s, 168 s, and 34 s for the 
IRAC bands and the MIPS 24, 70, and 160/im arrays, respectively. The data were reduced using 
the standard SSC pipeline and were coadded and corrected offline as needed. The details of the 
data reduction will be described in future FLS data papers (IRAC: Lacy et al. 2004; MIPS-24;um: 
Fadda et al. 2004b; MIPS-70,160;um: Frayer et al. 2004b). 



3. Results and Discussion 

3.1. Source Identification 

We observed 28 Westerbork radio sources with SCUBA and detected 7 SMGs at S/N> 3. 
Fourteen sources were not detected at 850^m, and the remaining 7 sources have marginal results, 
showing positive signals of 1.5-3o" (Column [4], Table 1). The Westerbork radio data have a resolu- 
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tion of 14" X 11" which is well matched to the SCUBA data, but is not sufficient to obtain reliable 
counterparts in general for the optical and Spitzer data. We used the higher resolution (5"0) VLA 
(1.4 GHz) data of the field (Condon et al. 2003) to obtain more accurate radio positions. The res- 
olution of the VLA data is well matched to the MIPS-24 resolution of 6". Although the VLA data 
have lower S/N (rms~ 20/iJy), the Westerbork sources were typically detected at > 3a in the VLA 
image, providing radio positional errors of about 1-2". The offsets between radio and MIPS-24 
positions are typically less than 2", consistent within the positional uncertainties of the radio data 
and the 1.5-2" positional uncertainties of the MIPS-24 data set. Based on the 24;^m counts in the 
FLS field (Marleau et al. 2004), the probability of a chance coincidence within 2" is only about 
1% for 24/im sources brighter than 100/uJy. Hence, confusion at 24/Lim is not a significant issue for 
this study. After the identification of the MIPS-24 and R-band counterparts based on the radio 
positions, the corresponding MIPS-24 and optical positions were used to identify the appropriate 
IRAC counterparts. 

Four of the Westerbork sources have multiple radio components at the resolution of the VLA 
data (1, 45, 79, 128). Source 1 is the brightest Spitzer source and is comprised of five components, 
all detected in the MIPS-24, R-, and the IRAC-bands. The details of the source identification are 
provided in Table 1. Figure 1 shows the MIPS-24 and IRAC-3.6 images for the 7 SMGs detected 
at 850/im with S/N> 3. 

All 14 sources showing 850;um emission above the 1.5a level have 24/im counterparts, although 
source 208 is only detected at about the 2.5-3(j level at 24^m. The detection of 24/im counterparts 
associated even with the marginal 850;um sources gives credence to the reliability of the > 3a 
SMG sources. In comparison for the 14 radio sources without SCUBA detections, only half have 
24^m counterparts. The detection of 24;um counterparts for this sample of radio-selected SMGs is 
consistent with the high fraction of 24^m counterparts found for the MAMBO and SCUBA sources 
in the Lockman Hole (Ivison et al. 2004; Egami et al. 2004). 

Seven sources in the full sample do not have optical counterparts down to R= 26.4 mag, and 
only three are not detected by IRAC. The only source without any optical or IRAC counterpart is 
source 108 which is the brightest radio galaxy in the sample (likely a radio-loud AGN). Two sources 
(199 [z = 1.06] & 150 [z = 0.84]) in the sample have spectroscopic redshifts from the ongoing Keck 
DEIMOS redshift survey (Choi et al. 2004). The low redshift of z = 1.06 for the 850^m source 199 
implies a cool dust temperature of 20 K assuming a temperature dependent sub-mm to radio redshift 
relationship (Blain 1999). 



3.2. Infrared Properties 

The MIPS 24^m flux densities for the galaxies detected at 850;um are consistent with expec- 
tations, assuming standard ULIRG spectral energy distributions (SEDs) at the typical redshifts of 
z ~ 2-3 found for the SMG population (Chapman et al. 2003b; 2004). Figure 2 shows a range 



of S24/S(1.4GHz) flux density ratios for the population that may reflect a wide distribution of 
infrared colors. All of the galaxies showing 850//m emission lie within the range of values expected 
for local ULIRGs redshifted to z = 1-4. Galaxies without 850/im emission show a wider range 
of S24/S(1.4GHz) flux density ratios. About half of the radio sources without detectable 850/im 
emission have properties consistent with ULIRGs, while the other half show excess radio emission 
compared to their 24//m emission, consistent with different degrees of radio-loudness. 

MIPS 70fiiQ observations can constrain the infrared colors of the SMGs. Infrared- warm, AGN 
dominated sources are expected to have flux density ratios of S70/S24 ~ 5, while starburst dom- 
inated, infrared-cool sources are expected to have ratios of S70/S24 > 10. We only detect two 
sources at 70/um, both of which are cool in the infrared (source 1: S70/S24= 11 and source 47: 
S70/S24= 26). To derive an estimated average S70/S24 ratio for the SMG population, we coad- 
ded the 70/im data at the radio positions for the sources showing 850/im emission. We find an 
upper-limit of S70< 1.2 mJy {2a). This corresponds to an infrared color of S70/S24< 5, leaving 
out the two sources with 70/im detections. Including the two 70/im detections, the average ratio 
for the SMG population is S70/S24 < 7, which is slightly lower than expected if the population is 
dominated by star-formation. Taken at face-value, these results could suggest that many SMGs are 
infrared-warm AGN, contrary to previous conclusions that the population is dominated by star- 
formation. However, the lower than expected S70/S24 ratios may arise from strong 7.7/im PAH 
emission redshifted into the 24/im band, given that the median redshift for the SMG population is 
z ~ 2.4 (Chapman et al. 2003b; 2004). 

The SEDs in the IRAC bands can be used to estimate redshifts from the rest-frame 1.6/im peak 
expected for star-forming systems (e.g., Egami et al. 2004) and to help identify AGN showing hot 
dust implied by their red IRAC colors (e.g., Ivison et al. 2004). For sources showing a significant 
bump in the IRAC bands, the implied redshifts are generally consistent (Column [11]&:[12], Table 
1) with the redshift estimates derived from the sub-mm to radio spectral index (Carilli & Yun 
1999; Blain 1999). The discrepant photometric redshifts for sources 48&: 119 may indicate cool dust 
temperatures (Td < 40 K) for these SMGs. 

Of the sources with 850/im emission and IRAC detections, only source 47 shows a strong 
increase in flux density as a function of wavelength across the IRAC bands consistent with an AGN. 
However, source 47 also has the highest S70/S24 ratio in the sample (S70/S24=26), suggesting that 
it is likely an infrared-cool starburst that could be at z > 4 and/or be highly reddened due to dust 
extinction. Hence, IRAC colors are not always conclusive for determining the properties at longer 
infrared wavelengths. 



3.3. Extremely Red Objects 

A significant fraction of SMGs is comprised of extremely red objects (EROs) (Small et al. 
1999; Frayer et al. 2004a). The ERO definition of i? — i^ > 6 (Thompson et al. 1999) corresponds 



to a flux density ratio of S3.6/S(R)> 85. Using this criterion, we identify 8 EROs (Fig. 3). The 
30% (8/27) fraction of EROs in this sample is larger than the 8% found for the total micro-Jy 
radio population (Small et al. 2002), presumably since the SCUBA targets in this sample were 
generally selected on the basis of being faint and red in the optical bands. Figure 3 shows a wide 
range of IRAC colors for the SMG population, and no correlation is observed between the IRAC 
colors and the R — [3.6] colors. Another interesting result is that most of the radio-selected EROs 
without detectable 850/um emission (6/7) are bright at 24^m, indicating that these sources are 
infrared-bright galaxies (likely high-redshift ULIRGs) below the current SCUBA detection limits. 



4. Conclusions 

All of the radio-selected SMGs in the sample have Spiteer counterparts, showing a wide range of 
infrared colors consistent with ULIRGs at z = 1-4. The combination of accurate radio positions and 
24^m detections is a powerful tool for the identification of SMGs. More sensitive observations are 
required in the MIPS-70 and MIPS-160 bands to measure the infrared colors of the SMG population. 
In addition, observations with the Spitzer InfraRed Spectrograph are needed to determine the 
level at which PAH emission contributes to the MIPS-24 flux densities and to help determine the 
AGN fraction of the population. In general, Spitzer selected ULIRGs will be biased toward AGN, 
except within specific redshift bins associated with the PAH features which may be biased toward 
starbursts. 

We thank the staff at the JCMT and the Spitzer Science Center for their support of these 
observations. This work is based in part on observations made with the Spitzer Space Telescope, 
which is operated by the Jet Propulsion Laboratory, California Institute of Technology under NASA 
contract 1407. 
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Fig. 1.— IRAC-3.6/im and MIPS-24/im images of the seven SMGs detected at 850//m with S/N> 3. 
Each panel is approximately 30" in size (north is up and east is left) with tickmarks every 1" . The 
grey-scale is plotted on a logarithmic scale, and the contours start at 3cj and increase by factors of 
2. The crosses represent the positions of the radio counterparts. 
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Fig. 2.— The S24/S(1.4 GHz) flux density ratios as a function of the S850/S(1.4GHz) flux density 
ratio. The solid circles are the > 3a SMGs, the open circles are the radio sources without 850/im 
emission, and the solid triangles are the sources with 850/im measurements of 1.5~3a. All limits are 
2a. The dashed lines show the expected ratios assuming power-law approximations for the infrared, 
sub-mm, and radio emission. The top dashed line represents an infrared-warm ULIRG, while the 
lower dashed line represents an infrared-cool ULIRG. Redshift estimates are shown by the dotted 
lines. The observed ratios are consistent with the SEDs of local ULIRGs at the expected redshifts 
of z --^ 1-4 for the SMG population. Strong radio galaxies which do not obey the far-infrared to 
radio correlation are located in the lower-left. 
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Fig. 3. — The S3.6/S(R-band) flux density ratios as a function of IRAC color given by the S8/S3.6 
flux density ratio for sources with detections at 3.6;uni. Symbols are the same as Figure 1. EROs 
are located above the dotted line, corresponding to i? — [3.6] > 7.5 (or R — K > 6). The data show 
a wide range of IRAC colors. 



Table 1. FLSV Radio-Selected Candidate Submillimeter Galaxies 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


Source 


Position 


S(1.4GHz) 


S(850Mm) 


S24 


_R-band 


S3. 6 


S4.5 


S5.8 


S8.0 


submm 


IRAC 


Notes 




a(J2000) 5(J2000) 


(MJy) 


(mjy) 


(MJy) 


(mag) 


(MJy) 


(MJy) 


(MJy) 


(/^Jy) 


(l + Zphot) 

(Td/40K) 


(1 + ^phot) 




1 


17:18:12.9+59:39:22 


750 ± 120 


7.0 + 2.3 


2750 + 280 


19.68 + 0.05 


559 + 56 


365 + 37 


362 + 36 


293 ± 29 


2.3 + 0.8 


< 2.5 


a 


44 


17:17:29.7+59:54:29 


74 ±9 


9.9 + 2.9 


399 + 40 


23.44 + 0.10 










4.4 + 0.8 




b 


48 


17:17:33.7+59:53:56 


54 ±9 


13.1 + 2.9 


232 + 28 


23.83 + 0.11 


16.5 + 1.7 


18.0 + 1.8 


14.9 + 3.3 


< 10 


5.0 + 0.8 


3 + 0.8 




115 


17:17:43.4+59:48:03 


262 ± 26 


5.4 + 1.6 


373 + 37 


24.72 + 0.20 


13.0 + 1.3 


13.0 + 1.3 


< 10.0 


28.1 + 3.3 


2.8 + 0.8 






119 


17:17:12.7+59:47:53 


81 ±9 


10.2 + 2.8 


146 + 28 


>26.4 


43.3 + 4.3 


44.0 + 4.4 


42.9 + 4.3 


25.5 + 3.3 


4.3 + 0.8 


3 + 0.8 


c 


199 


17:17:29.4+59:41:13 


116 ±12 


7.5 + 2.5 


301 + 30 


23.44 + 0.14 


35.4 + 3.5 


39.9 + 4.0 


49.3 + 4.9 


48.6 + 4.9 


3.7 + 0.8 




z = i.oe'' 


208 


17:18:10.9+59:40:41 


77 ±9 


6.8 + 1.9 


78 + 28 


22.62 + 0.10 


24.5 + 2.5 


24.0 + 2.4 


< 10 


19.3 + 3.3 


3.9 + 0.8 






45A 


17:17:47.5+59:54:24 


80 + 9 


7.2 + 2.5 


200 + 28 


24.12 + 0.14 


26.0 + 2.6 


27.4 + 2.7 


36.5 + 3.7 


47.4 + 4.7 


4.0 + 0.8 




c 


47 


17:17:22.5+59:54:12 


100 ± 10 


4.9 + 2.4 


498 + 50 


25.31 + 0.22 


28.2 + 2.8 


45.6 + 4.6 


68.1 + 6.8 


111.6 + 11.2 


3.4 + 0.8 




c,f 


73 


17:17:57.9+59:52:00 


118 + 12 


5.8 + 2.1 


135 + 28 


> 26.4 


29.4 + 2.9 


35.5 + 3.6 


30.7 + 3.3 


27.2 + 3.3 


3.4 + 0.8 


3 + 0.8 


c 


139 


17:18:23.2+59:45:53 


215 ± 22 


5.0 + 2.8 


444 + 44 


22.62 + 0.05 


36.9 + 3.7 


32.5 + 3.3 


27.9 + 3.3 


24.1 + 3.3 


2.9 + 0.8 


<2.5 




145 


17:17:46.2+59:45:17 


66 ±9 


8.7 + 4.0 


161 + 28 


>26.4 


7.2 + 1.0 


7.7+1.0 


< 10 


19.1 + 3.3 


4.3 + 0.8 




c 


156 


17:18:16.8+59:44:30 


82 ±9 


4.4 + 2.5 


368 + 37 


25.64 + 0.37 


9.7+1.0 


12.9 + 1.3 


25.0 + 3.3 


21.0 + 3.3 


3.5 + 0.8 


3.5 + 0.8 




191 


17:17:15.5+59:42:02 


145 ± 15 


3.8 + 2.5 


277 + 28 


25.35 + 0.34 


<3 


<3 


< 10 


< 10 


3.0 + 0.8 






75 


17:18:01.7+59:51:47 


38 ±9 


<7.5 


416 + 42 


> 26.4 


20.0 + 2.0 


24.6 + 2.5 


27.2 + 3.3 


19.8 + 3.3 


<4.8 


3.5 + 0.8 


c 


79A 


17:17:22.8+59:51:30 


114 + 11 


<6.6 


408 + 41 


23.30 + 0.17 


54.2 + 5.4 


38.9 + 3.9 


36.6 + 3.7 


25.7 + 3.3 


<3.6 


< 2.5 


g 


85 


17:18:12.3+59:50:56 


44 ±9 


<8.4 


<84 


21.93 + 0.05 


40.1 + 4.0 


26.4 + 2.6 


23.9 + 3.3 


16.5 + 3.3 


<4.7 


<2.5 




91 


17:17:43.7+59:50:22 


43 ±9 


<7.8 


331 + 33 


> 26.4 


13.3 + 1.3 


17.9 + 1.8 


15.5 + 3.3 


18.3 + 3.3 


<4.7 




c 


99 


17:17:06.4+59:49:25 


139 ± 14 


<8.4 


< 84 


23.50 + 0.15 


36.3 + 3.6 


25.3 + 2.5 


10.0 + 3.3 


< 10 


<3.6 


<2.5 


h 


108 


17:17:41.2+59:48:36 


1080 ± 130 


< 8.1 


<84 


> 26.4 


<3 


< 3 


< 10 


< 10 


<2.2 






109 


17:17:38.5+59:48:32 


765 ± 110 


< 6.6 


<84 


24.19 + 0.14 


<3 


< 3 


< 10 


< 10 


< 2.2 






128B 


17:17:47.0+59:47:12 


62 ±9 


<9.0 


< 84 


25.41 + 0.26 


4.1 + 1.0 


3.6 + 1.0 


< 10 


< 10 


< 4.4 




i 


155 


17:17:56.0+59:44:32 


60 + 9 


< 7.2 


<84 


24.60 + 0.27 


30.0 + 3 


32.7 + 3.3 


21.7 + 3.3 


< 10 


<4.3 


3 + 0.8 




136 


17:17:11.6+59:46:21 


108 + 9 


<6.9 


292 + 29 


23.89 + 0.16 


21.8 + 2.2 


23.5 + 2.4 


19.0 + 3.3 


17.3 + 3.3 


<3.7 


3 + 0.8 




146 


17:18:12.9+59:44:54 


627 ± 63 


<9.9 


217 + 28 


23.83 + 0.10 


39.3 + 3.9 


27.4 + 2.7 


15.7 + 3.3 


< 10 


< 2.6 


< 2.5 




150 


17:17:42.4+59:44:56 


68 ±9 


<6.9 


259 + 28 


23.37 + 0.14 


18.7 + 1.9 


11.8 + 1.2 


< 10 


21.8 + 3.3 


<4.1 




z = 0.84'^ 


198 


17:17:21.4+59:41:13 


114 + 12 


<8.1 


< 84 


> 26.4 


21.2 + 2.1 


22.0 + 2.2 


13.3 + 3.3 


13.8 + 3.3 


<3.7 


< 2.5 


c 


211 


17:18:21.3+59:40:27 


356 ± 36 


<8.4 


432 + 43 


26.16 + 0.42 


32.3 + 3.2 


41.7 + 4.2 


60.1 + 6.0 


89.2 + 8.9 


< 2.9 




c 



Note. — (l)Radio source identifier. (2)Radio positions which are accurate to about +2". (3)Radio flux densities from the Westerbork and VLA data. (4)SCUBA 
850/im measurements. Sources 1—208 have > 3cr 850/im detections. Sources 45A-191 have 1.5— 3<t 850/im measurements, and sources 75—211 are non-detections at 
850/im. All flux density and magnitude hmits are given as 3it. (5)MIPS-24/im measurements. (6)0ptical ij-band magnitudes from NOAO 4m. (7)— (lO)IRAC flux 
densities. (ll)Photometric redshift estimate of (1 + z) X {T^/40K)~^ based on the sub-mm to radio relationship of CariUi & Yun (1999), accounting for the degeneracy 



of redshift and dust temperature in the relationship (Blain ct al. 1999). (12)Photometric redshift based on the assumed rest-frame 1.6^m peak in the IRAC SED. 
(a)Measurenients summed over all five components (Fig. 1). Source 1 is also detected at longer wavelengths with S70= 30 it 6mjy and S160= 130 it 26mjy. Sources 
not detected in the 70&:160/im MlPS-bands have 3a limits of about < 10&< 50mjy, respectively. (b)Detected in IRAC-bands, but IRAC measurements are confused 
with nearby bright optical source. (c)ERO based on high S3.6/S(R) flux density ratio, (d) Spectroscopic redshift based on Keck DEIMOS data (Choi et al. 2004). 
(e)Brighter radio component 45A has 24/im counterpart and is adopted as the primary counterpart, since 45B has no associated 24fiin source. (f)Detected in MIPS- 70 
with S70= 13 it 3mjy. (g)Brighter radio component 79A has a 24/im counterpart and 'merger-like" morphology with two optical components. (h)Sourcc99 has two 
optical components showing a merger-like morphology. (i)Neither radio component of source 128 has a 24^m counterpart; 128B is adopted as the counterpart for flux 
density comparisons since there are no optical or IRAC counterparts associated with 128A. 



I 



